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We investigate the origin of superconductivity in doped SrTiC >3 (STO) using a combination of 
density functional and strong coupling theories within the framework of quantum criticality. Our 
density functional calculations of the ferroelectric soft mode frequency as a function of doping reveal 
a crossover from quantum paraelectric to ferroelectric behavior at a doping level coincident with 
the experimentally observed top of the superconducting dome. Based on this finding, we explore 
a model in which the superconductivity in STO is enabled by its proximity to the ferroelectric 
quantum critical point and the soft mode fluctuations provide the pairing interaction on introduction 
of carriers. Within our model, the low doping limit of the superconducting dome is explained by the 
emergence of the Fermi surface, and the high doping limit by departure from the quantum critical 
regime. We predict that the highest critical temperature will increase and shift to lower carrier 
doping with increasing 18 0 isotope substitution, a scenario that is experimentally verifiable. 


Strontium titanate (STO) is a cubic perovskite with 
the ideal prototype structure at room temperature and 
a tetragonal structure below ~100K due to symmetry¬ 
lowering antiferrodistortive (AFD) rotations of the TiOg 
octahedra. It is characterized by a number of remarkable 
properties. It was the first superconducting oxide to be 
discovered and shows a dome as a function of doping, sim¬ 
ilar to that of the high-T c cupratesp], but with its max¬ 
imum transition temperature at T c ~ 0.4K. Early tun¬ 
neling measurements [2] and subsequent experiments [3] 
suggested an unusual two-band superconductivity, con¬ 
sistent with the closely spaced lowest conduction bands 
at the center of the Brillouin zone. In addition, the onset 
of superconductivity has been shown to occur at remark¬ 
ably low carrier concentrations of 10 18 e/cm 3 [3j. De¬ 
spite a long-running interest in its origin T], a complete 
theoretical account of the superconducting dome remains 
elusive, and many aspects of superconductivity in STO 
remain a puzzle. 

The dielectric behavior of STO is also unusual. The 
dielectric constant is strongly temperature dependent, 
and diverges at low temperature in a manner character¬ 
istic of a ferroelectric phase transition [2j. Rather than 
manifesting ferroelectric behavior, however, STO is a so- 
called quantum paraelectric, in which quantum fluctua¬ 
tions at zero temperature suppress the transition to the 
ferroelectric state [3]. The quantum paraelectric state is 
characterized by low energy excitations and large ferro¬ 
electric fluctuations |5], and it has been speculated that 
these might be relevant for the superconductivity 
Indeed, early descriptions mm of the superconducting 
dome in STO were based on the effects of screening of the 
interaction between electrons and the optical phonons re¬ 
sponsible for the large dielectric response. Because heav¬ 
ier 18 0 atoms suppress the quantum fluctuations, STO 
develops ferroelectric order on isotope substitution of 16 0 
with 18 0, and the composition with 35% 18 0 substitu¬ 


tion was recently reported to be a ferroelectric quantum 
critical point (QCP) [9]. 

We present a model in which these two features - prox¬ 
imity to the ferroelectric QCP and the unusual supercon¬ 
ducting properties - are intimately related, and the su¬ 
perconducting dome emerges as a result of the quantum 
critical ferroelectric fluctuations. A connection between 
the formation of a superconducting dome and quantum 
criticality has been extensively discussed in the context of 
unconventional superconductivity, both in heavy fermion 
materials and in the cuprates EIDHH]- It is proposed that 
competing phases close to the quantum critical point lead 
to low energy excitations such that any residual interac¬ 
tions drive the system to a new, possibly superconduct¬ 
ing phase. In heavy fermion materials and the cuprates a 
magnetic quantum critical point with associated spin ex¬ 
citations has been invoked to explain superconductivity. 
In STO the elementary excitations associated with the 
ferroelectric quantum critical point are optical phonon 
modes. As a result we expect differences in the na¬ 
ture of the superconducting order: Magnetic fluctuations 
typically produce unconventional superconducting order 
such as d -wave singlets for antiferromagnetic fluctuations 
mm or p- wave triplet states for ferromagnetic fluc¬ 
tuations m ■ The ferroelectric fluctuations in STO, in 
contrast, involve q = 0 phonon modes and as such are 
candidates for pairing interactions that introduce con¬ 
ventional s-wave superconducting order, as observed in 
experiments }16| . 

Our model is motivated by our density functional the¬ 
ory (DFT) calculations of the zone-center (q = 0) soft- 
mode optical phonon frequency as a function of electron 
doping, which shows an intriguing correlation with ex¬ 
perimental measurements of the superconducting dome. 
This ferroelectric soft mode, which consists predomi¬ 
nantly of opposite Ti cation and O anion displacements 
(for details, see Ref. HTJ), has a calculated imaginary fre- 
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FIG. 1. Literature values of the superconducting critical tem¬ 
perature [I] (circles) and calculated frequencies (this work) of 
the ferroelectric modes parallel- (||) and perpendicular (_L) to 
the axis of the AFD rotations (red solid and dashed lines) as 
a function of the carrier concentration. The imaginary fre¬ 
quencies obtained at low doping indicate negative restoring 
forces corresponding to ferroelectric instabilities; as the car¬ 
rier concentration is increased the ferroelectric mode hardens 
and its phonon frequency becomes real. The inset shows the 
calculated energy as a function of ferroelectric mode ampli¬ 
tude for various doping levels, illustrating the crossover from 
the classic ferroelectric double well potential energy to a sin¬ 
gle well, indicating a paraelectric ground state on increasing 
doping. As the charge carrier concentration is increased, T c 
first increases and then decreases, forming the characteristic 
superconducting dome. We see that the doping concentration 
at which T c drops to zero, ~ 10 2O e/cm 3 , closely matches that 
at which the ferroelectric mode hardens. 


quency at zero doping, indicating the presence of a fer¬ 
roelectric instability. The calculated potential energy as 
a function of the relative position of anions and cations 
(see supplementary material Sec. I for details) shows 
the characteristic double well form, with the two minima 
corresponding to ferroelectric structural ground states 
with opposite polarizations. In practice, quantum fluctu¬ 
ations between the two wells suppress the ferroelectricity 
in STO, and give it its quantum paraelectric behavior. 
In Fig. [T] we show how, on electron doping, the mod¬ 
ulus of the mode frequency decreases, corresponding to 
a weakening of the ferroelectric instability, and the fre¬ 
quency eventually becomes real - signaling a single high- 
symmetry energy minimum (inset to Fig. [l]) - at a dop¬ 
ing concentration of ~ 10 2 ° cm' 3 . Since there is now 
only one minimum of the potential well, there are clearly 
no quantum fluctuations between equivalent states. At 


the same doping level, the experimentally measured su¬ 
perconducting transition temperature starts to reduce. 
Since soft modes are characteristic of quantum critical¬ 
ity [T§] , we propose therefore the following model for the 
superconducting dome in STO: First, superconductivity 
is favored when the quantum fluctuations favored by the 
soft lattice modes increase the superconducting coupling 
constant A. However, these are strongest at low doping, 
where there are insufficient carriers to provide robust su¬ 
perconductivity. Increasing the doping level has the side- 
effect of reducing A, which in turn determines the upper 
bound of the superconducting dome. 

To test this hypothesis we propose isotopic substitu¬ 
tion of 16 O with 18 O, which lowers the energies of the 
zero-point energy levels in the two minima and reduces 
the probability of tunneling between them eventually fa¬ 
voring a ferroelectric ground state, see Fig. n>. We know 
that at zero carrier doping, the paraelectric to ferroelec¬ 
tric transition occurs at about 35% 18 O substitution and 
is a quantum critical point [9]. In addition, our DFT cal¬ 
culations tell us that doping suppresses FE, and so the 
QCP should move to higher 18 O fractions as doping is 
increased, implying the existence of a quantum critical 
line (QCL). This allows us to construct the schematic 
phase diagram in Fig. [2j Our DFT calculations give an 
upper bound for this QCL, which is the doping level at 
which the frequency of the FE mode becomes real and 
the quantum fluctuations are completely suppressed; in 
practice this represents the limit of infinitely heavy oxy¬ 
gen atoms and the actual critical transition will occur 
at much lower doping. Note that, at least at low 18 0 
concentration, charge carriers only appear as one moves 
away from the quantum critical point, so the QCP is in 
fact located outside the superconducting dome. This is in 
contrast to the emergence of superconductivity in other 
systems close to a QCP, such as the cuprate supercon¬ 
ductors, in which the dome is approximately centered on 
the QCP. In those cases, the QCP occurs at substantial 
doping, where charge carriers are already available. 

We quantify our proposed model by calculating T c , as¬ 
suming the scenario of soft critical modes in the limit 
of low doping. We first write a quantum model for the 
ferroelectric phase transition which yields a spectrum for 
the FE phonons. Then, we calculate the superconducting 
coupling constant, using the Macmillan formula [2D] . 

We use the order-disorder approach m to model the 
ferroelectric fluctuations of the modes shown in Fig. [T] 
We assume that these modes have Ising character. By 
analogy with magnetic phase transitions, the transverse 
Ising model 


i i,j 

can be used to describe the FE transition [5]. Here cr x ,z(i) 
are the Pauli matrices for site i, T/h is the onsite tunnel- 
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FIG. 2. a) Schematic phase diagram of STO as a function 
of carrier doping and isotope replacement. The orange cir¬ 
cles mark the experimentally measured transition to super¬ 
conductivity, as observed in Ref. [0 The blue circles are the 
measured transition temperatures m from the paraelectric 
(PE) to the ferroelectric (FE) phase as a function of 18 O iso¬ 
tope substitution. Our DFT calculations suggest that the 
ferroelectric phase penetrates slightly into the non-zero dop¬ 
ing regime, but then quickly disappears as doping suppresses 
ferroelectricity, although no experimental data for this tran¬ 
sition line is available. The maximal value of doping at which 
the ferroelectric phase persists is labeled as n*. Although we 
have no precise calculation for n* , its value should lie in the 
range 10 19 < n* < 10 2l] . b) Schematic illustration for the 
lowering of the lowest energy levels (dashed red lines) in the 
double well potential (black solid line) as /is is increased. 


ing rate, Jij is the inter-site coupling, given by the energy 
difference between two cells with their dipoles aligned 
parallel or anti-parallel to each other, and the eigenstates 
of a z represent the state of the system in one of the 
two wells. The quantum phase transition occurs when 
T ~ Jo,j [5]. Our DFT study shows that doping the 
system will reduce the barrier and thus increase T. The 
excitations of |l|, in the paraelectric phase T > J2j Jo,j> 
are given by |5j 


w2 = 4r(r-<«7 X >j q ) (2) 

where J q = ]T) J 0 j e* Rjq is the Fourier transform of the 
coupling and (<j x ) ~ 1 is the average of a x (i). In our 
analysis we consider only nearest-neighbor coupling for 
simplicity. Long range interactions make the calculation 
more intricate but do not yield any qualitative changes. 
Furthermore, since the antiferrodistortive rotations of the 
TiOe octahedra render the lattice highly anisotropic, we 
treat the system as one dimensional. Thus, we write the 


coupling as J q = 2 J cos(g), where J is a constant and q is 
the wave number in the direction of the largest coupling. 

When the system is close to the phase transition it be¬ 
comes gapless as the lowest excitation softens, w q=0 —> 0 
(see supplementary material, Sec. II)). This is accom¬ 
panied a large susceptibility and an enhanced electron- 
phonon coupling. To quantify this idea we calculate the 
dependence of T c on the phononic spectrum using the 
formalism of Eliashberg strong-coupling theory. The cou¬ 
pling constant for superconductivity is given by [20] 

A = f a 2 (uj)F{uj) — (3) 

Jo w 

where a(cu) is the electron-phonon coupling, which we 
assume to be the constant a, and F(ui) is the spectral 
density of the phonons. In the limit of a van Hove singu¬ 
larity at q = 0, so that F(u) ~ S(ui — oj o), this yields 


w q =o(/i8> Ef) 


which already captures the main physical picture of soft- 
mode enhanced superconductivity. The full solution is 
obtained by inserting F(ui) = f dqS(ui — ui q ) into ^ and 
transforming it to an integral over q: A = f a 2 where 
uiq is given by <©• One then obtains 


a ~ r dq 

J-tt 2T— cos( 9 )/r 


(5) 


The critical temperature can then be obtained by com¬ 
bining this coupling constant with the standard expres¬ 
sion (see for example Ref. [2^1) 




deN(e) 


tanh (e/2 T c ) 
e 


( 6 ) 


where e is the energy relative to the Fermi energy, Ep, 
and N(e) is the density of states. The lower limit of the 
integral is set by N(e) = 0 at and below the bottom of 
the band where e < —Ep. The upper limit is set by 
the Fermi level, where we define e = 0. Since in the low 
doping scenario that we consider here the relevant energy 
range is close to the bottom of the band, we can assume 
that N(e) ~ y/e + Ep close to e = — Ep. Using x = e/T c 
equation § then becomes 


D 

X 


vfj" < 7 ) 

J-E f /T c x 


where D is a constant of proportionality. Note that T c 
has a double dependence on Ep: one directly from the 
limit of the integral in (J7| and the other from the de¬ 
pendence of A on the tunneling rate T on Ep through its 
dependence on the carrier concentration. 

Before we can solve Eq. [T] numerically to obtain T c as 
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FIG. 3. Calculated T c as a function of doping level for several 
fractions of isotope replacement, fig. The blue diamonds are 
experimental results taken from Ref. JT|. Replacing 16 O with 
18 0 moves the QCP closer to the doping range relevant for 
superconductivity and causes a significant enhancement in T c . 
we use the parameters A = 0.4, B = 10~ 6 K -2 , C = 2.5 x 
10~ 3 K _1 , D = 95 K 1,/2 , as defined in the main text. 

a function of Ep , we need the explicit dependence of the 
parameters of our model on doping and isotope re¬ 
placement. The quantity with the largest quantitative 
influence is the ratio T /2 J, which is equal to one on the 
QCL. For simplicity we set 2J = 1 and consider only 
the dependence of T. As discussed above, carrier dop¬ 
ing decreases the barrier between the two wells and thus 
increases the tunneling energy F. This effect starts at 
low doping level and becomes very strong around carrier 
concentrations of 10 2O e/cm 3 , so we consider both a lin¬ 
ear and quadratic dependence on Ep- 18 0 replacement 
on the other hand should decrease T approximately lin¬ 
early as the zero-point energy levels shift deeper into the 
wells. Furthermore, we require that at zero doping and 
35% 18 0 substitution, which is the known QCP, F should 
equal unity. Thus we use the form 

r = 1 - A(f ls - 0.35) + BE 2 p + CEp, (8) 

where / 18 is the 18 0 fraction and the constants A, B, 
and C are chosen so that the calculated T c for /i 8 = 0 
matches the experimental value. We then use the ex¬ 
pression r from Eq. ([8]) and insert this into Eq. (J7| to 
calculate T c . 

In Fig. [3] we plot our calculated T c as a function of the 
Fermi energy (converted to carrier concentration) for var¬ 
ious values of fig. Two features are clear from the plot: 
i), we find a significant enhancement of T c with increased 
18 O content, reflecting the fact that the isotope substi¬ 
tuted system is closer to the QCP. ii), we find that the 
peak of the superconducting dome shifts to lower carrier 
concentrations, since the enhancement of A and thus T c is 
strongest close to the QCP, as can be seen from Eq. (|3]). 
We note that, even when fig exceeds 0.35, doping quickly 
reduces the depth of the double wells, allowing quantum 


fluctuations to return STO to the quantum paraelectric 
state. Thus, apart from the limit of very low doping, all 
systems we consider have paraelectric, not ferroelectric 
ground states. In our mechanism for superconductivity 
in STO, increasing the atomic mass leads to an increase 
of the critical temperature. That is /T c > 0 (for de¬ 
tails, see supplementary material, Sec. III). This differs 
profoundly from the well-known isotope effect in BCS su¬ 
perconductors, in which P2j, where M is 

the mass of the atoms. This arises from the dependence 
of T c on the Debye frequency. 

We have provided a description of the superconduct¬ 
ing dome in STO in which the QCP at zero doping pro¬ 
vides low energy soft phonon excitations, which lead to 
a large coupling constant. Increasing the doping pro¬ 
vides carriers for superconductivity but reduces the fer¬ 
roelectric quantum fluctuations and decreases the cou¬ 
pling constant, eventually suppressing the superconduc¬ 
tivity and limiting the top of the superconducting dome. 
Since isotope substitution allows tuning of the QCP, our 
model predicts a large and unusual isotope effect on T c , 
see Fig. [3j which should be experimentally observable. 
The understanding of the competition between carrier 
concentration and proximity to a QCP developed here 
provides a new design guideline in the search for novel 
superconducting compounds and suggests a route to en¬ 
gineering materials with higher T c s through tuning the 
location of their QCP. 
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FIG. 4. The displacements of the atoms within the unit cell for a) the AFD mode as well as b) the FE mode parallel to the 
AFD axis and c) the FE mode perpendicular to the AFD axis. 
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Supplemental material 
Ferroelectric modes 

The ferroelectric mode frequencies plotted in Fig. 1 in the main text refer to the mode parallel to the axis of the 
antiferrodistortive rotations of the Ti0 6 octahedra (solid line) and to the doubly degenerate perpendicular modes 
perpendicular to the octahedral rotations (dashed line). The displacements of the atoms within the unit cell for each 
of these modes are shown in Fig. [4j 

The mode frequencies were computed using the frozen phonon approach as implemented in Phonopy [Tj with 
forces obtained from density functional theory with the PBEsol functional [2J as implemented in the Vienna ah initio 
simulation package (VASP) [3HB]. In these force calculations, wavefunctions were expanded in plane waves up to a 
kinetic energy cutoff of 550 eV for PAW potentials El with Sr(4s, 4p, 5s), Ti(3p, 3d, 4s) and 0(2s, 2p) states 
in the valence. Reciprocal space was sampled using a 4x4x4 Monkhorst-Pack mesh [9] for the 2x2x2 supercell used 
to describe the structure containing the fully relaxed AFD. Electron doping was performed by adding extra charge 
carriers to the cell while applying a neutralizing background charge. 


The frequency of the soft modes 

Our model predicts some unusual effect with regards to isotope replacement. One of them is the softening of the 
lowest frequency mode. This frequency, u; q= o, is given by Eq. (2) in the main text. By using our model for the 
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FIG. 5. The fractional change in the superconducting critical temperature, /jU , as a function of 18 0 fraction, for several 
carrier concentrations. The plots are based on the numerical calculations used to produce the superconducting dome in Fig 3 
in the main text 


dependence of the tunneling energy on the 18 0 fraction (Eq. (8) in the main text) we can explicitly calculate the 
change of this frequency with respect to the 18 O fraction (/is): 


du>q = o _ ^ 2r — 1 

dfia ~ ~ v^r- !)’ 


(9) 


For r > 1 this thus always lowers the mode frequency as the 18 0 fraction is increased, taking the system closer to the 
QCP. 


Equation for the change of T c 


Using Eq. (7) in the main text, we can calculate the change in T c as a function of /is- We obtain 


1 dT c 
T c d/ 18 


A , Act 1 2 fTT r 

2A + 4A 2 J-iT (r 2 -r COS q ) 3 4 / 2 
Ef tanh(ic/2) 1 

VTcD J-E f /T 0 2 xy/x+Ep/Tt 2X 


( 10 ) 


A careful study of this expression shows it to be positive for T > 1, as both the numerator and denominator are 
positive. More insight into this function can be gained from its plot, see Fig. [5] Apart from the fact that the change 
of T c is positive, we can see that it is larger for lower doping, i.e. when the system is closer to the QCP. This explains 
the shift of the peak of the dome to lower doping with increasing 18 0 fraction. For very small doping we can see that 
(dT c / dfis)/T c increases for large /is, implying a non-linear dependence on the mass. This suggests that the system 
is very close to the QCP, in contrast to the case of moderate and large doping, where the change is roughly constant 
for the available mass variation. 
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